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Abstract: A total of 44 taxa were observed in monthly sampled phytoplankton of the Shibpukur pool in Burdwan, West Bengal between 
March, 2010 and February, 2011. The most abundant taxa belong to Charophyta, followed by cyanobacteria, diatoms and euglenoids. 
Bio-indication shows that the pool community preferred low alkaline, low mineralized and low organically polluted water. The total 
phytoplankton density showed its maximum values in May, 2010. The Shannon -Weaver diversity index and the Pielou evenness value 
were found to be highest during the post-monsoon season. The Simpson dominance index and the Margalef index of richness were 
highest in the pre-monsoon season. The total phytoplankton density showed a highly significant positive correlation with pH and 
salinity and significant positive correlation with air temperature, water temperature, dissolved oxygen and total suspended solids. 
Among the reported 44 phytoplankton taxa only 3 showed a random distributional pattern. The Bray-Curtis Cluster analysis and the 
comparative statistics reveal two groups of phytoplankton assemblages in respect to the monsoon seasons. The successive communities 
form a continuum corresponding to Colwells Constancy (C) category. The calculated indices, CCA, and bio-indication analysis exhibit 
a low pollution level in the Shibpukur pool that can be used as a model of aquatic community dynamics under seasonal fluctuation in 
the monsoon climate, applicable for monitoring of water bodies in the West Bengal Province. 
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Introduction 

Temporal changes in the composition of phy- 
toplankton are of utmost importance to the aquatic 
systems health. The planktonic algae are primary pro- 
ducers and are at the base of the trophic cascade with 
zooplankton, shell fish, fin fish and with other higher 
aquatic animals. So a minute change in the composi- 
tion of the phytoplankton community can modify the 
food-web structure and the energy flow through the 
aquatic ecosystem. The phytoplankton contains di- 
verse prokaryotic and eukaryotic organisms. Accord- 
ing to Wetzel (1975) well-developed phytoplankton 
communities are normally restricted to stagnant wa- 
ter bodies and large slow rivers, whereas higher flow 
rates can disturb and disintegrate the phytoplankton 
structure. Various interacting environmental fac- 
tors regulate the temporal changes in the composi- 
tion of phytoplankton, its stability and predictability 



of changes (Colwell 1974). Temperature and light are 
climatic dependent variables, and are the key drivers 
which help to establish phytoplankton in the euphotic 
zone, as well as controlling availability of micro and 
macro nutrients and succession strategies (Reynolds 
1984; Harris 1987). Species of phytoplankton can be 
used as indicators of water quality (Palmer et al. 1977; 
Shubert 1984). They also play a decisive role in the 
global biogeochemical cycling (Broecker 1974). 

Though species composition and the seasonal 
dynamics of phytoplankton of stagnant water bod- 
ies of India have been investigated by several authors 
during the last century (Biswas 1949; Abraham 1962; 
Jana et al. 1980; Unni 1982; Verma et al. 1984; Zutshi 
1989; Gopal and Zutshi 1998; Ghosh and Keshri 2010, 
2011; Senapati et al. 2011), the relationships between 
algal communities' structure and environment is in- 
sufficiently studied (Ghosh et al. 2012) in respect to 
the monsoon climatic conditions. 
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The purpose of this study was to elucidate the 
changes in the composition of phytoplankton as a 
temporal variable of the influence of the monsoon cli- 
mate. Phytoplankton composition, its seasonal chang- 
es and the diversity characteristics of the Shibpukur 
pool have not been previously explored. The studied 
water body, typical for this area, is almost undisturbed 
by any anthropogenic sources and is partly covered by 
trees and herbaceous vegetation. The area is located 
adjacent to the Ramnabagan Wildlife Sanctuary, Burd- 
wan province, West Bengal. The environmental factors 
which are the principal contributors to the temporal 
fluctuation of phytoplankton communities are char- 
acterized by ecological indicators obtained in the 
course of the community structure analysis, involving 
various diversity indices and statistical approaches. 

Description of the study site 

In the West Bengal Province, the River Ganga 
valley abounds in small pools that occur in urbanized 
areas. The studied pool is located at 23° 15' 14.5" N, 
87°5ri0.5" E and 20 metres above sea level. This pool 
is rectangular in shape 50 m long and 25 m wide, shal- 
low (to one m deep) and has a very small catchment 
basin area because all the surroundings are urban- 
ized, except a 50 m zone around the pool itself (Fig. 1). 
The water body under consideration is located in the 
Burdwan district of West Bengal. It is locally known 
as "Shibpukur" of Shiva Temple Natural Reserve. Now 
only this water body and the temple are maintained 
by Burdwan Maharaja Trusty. The whole area is within 
a natural forest patch developed during the periods 
of the Rajas and Maharajas of Burdwan Taluk (Sarkar 
2008) and after the implementation of the Estate Ac- 
quisition Act, it was handed over to the West Bengal 
Forest Department for proper management and main- 
tenance. 

Experimental procedures 

To study the temporal changes in the phyto- 
plankton composition of a water body in the Burdwan 
district of West Bengal, India, monthly samplings from 
March 2010 to February 2011 were undertaken. Here 
for the construction of a seasonal framework we have 
designated the months March, April, May as summer; 
June, July, August as pre-monsoon; September, Octo- 
ber, November as monsoon and December, January, 
February as post-monsoon season. The geographic 



coordinates of the water body under study were de- 
termined by the GARMIN GPS map 76 CSx device. 
The location map was constructed with the help of 
ArcGis 10.0 software. 

In parallel with phytoplankton sampling, meas- 
urements were taken of air and water temperature by 
thermometer, electrical conductivity, total dissolved 
solids, salinity, dissolved oxygen and pH of water by 
Oakton waterproof Multiparameter Tester PCS 35 and 
Oakton portable DO Meter. 

Phytoplankton samples were collected in 400 
ml amber coloured bottles between 9 to 10 am and 
preserved with acetic Lugols solution (Sournia 1978). 
The collected samples were left for two days for sedi- 
mentation. Then the supernatant was withdrawn us- 
ing a pipette, leaving the settled cells in the bottom. 
The drop count method (Trivedy and Goel 1984) was 
practised for quantitative estimation of phytoplank- 
ton. Phytoplankton was studied with an Olympus GB 
microscope under magnifications of 960-1500 and 
for taxonomic identification prism drawings were 
done. Phytoplankton densities were calculated as cells 
per litre. The algal abundances were scored for indi- 
ces calculation (Table 1). For identification we used 
handbooks (Cox 1996; Desikachary 1959; Hustedt 
1930; Krammer and Lange-Bertalot 1991, 1997a,b,c; 
Komarek and Anagnostidis 1998, 2005; Prescott 1962; 
Smith 1950; Turner 1982; Wehr and Sheath 2003). The 
list of revealed species was arranged alphabetically ac- 
cording to M.D. Guiry and G.M. Guiry (201 1) updated 
taxonomic system. 



Table 1. Phytoplankton species abundance score scale (Barinova 
et al. 2006) 



Score Visual estimate Cell numbers per slide 



1 


Occasional 


1-5 


2 


Rare 


10-15 


3 


Common 


25-30 


4 


Frequent 


1 cell over a slide transect 


5 


Very frequent 


Several cells over a slide transect 


6 


Abundant 


One or more cells in each field of view 



The ecological and geographical characteris- 
tics of algal species were obtained from the database 
compiled for freshwater algae from multiple analyses 
of algal biodiversity (Barinova et al. 2006) with addi- 
tions (Ter Braak and Van Dam 1989; Van Dam et al. 
1994) according to substrate preference, temperature, 
streaming and oxygenation, pH, salinity, organic en- 
richments, N-uptake metabolism, and trophic states. 
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Fig. 1. Map of the study site of the Shibpukur pool in the Burdwan district, West Bengal 
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The ecological groups were separately assessed 
according to their significance for bio-indications. 
Species that respond predictably to environmental 
conditions were used as bio-indicators for particu- 
lar variables of aquatic ecosystems, the dynamics of 
which is related to environmental changes. Indices of 
saprobity and water quality classes were calculated af- 
ter Sladecek (1973). 

The specific diversity calculation was deci- 
phered from Shannon and Weaver (Shannon and 
Weaver 1949) as follows: 

where: H is the species diversity value, and p. is the 
proportion of individuals found in its species. 

Specific dominance was calculated from Simp- 
sons equation (Simpson 1949): 



^ = Z 



^7V[7V-1] 



J 



where: D is the specific dominance value, n is the 
number of individuals in its species, and N is the total 
number of individuals. 

The evenness index was calculated by the Shan- 
non evenness measurement process (Pielou 1969, 
1975) in the following way: 

J'=H'/H max =H'/\nS, 



where: / c is the evenness index, and H is the maxi- 

J max 

mum diversity that could possibly be found in a con- 
dition where all species had equal abundance values 
i.e. 



H'=H mm =lnS 

Species richness index was calculated following 
Margalef s diversity index (Clifford and Stephenson 
1975): 

D MG ={S-l)/lnN , 

where: D MG is Margalef s diversity index, S is the 
number of species found, and N is the total number of 
individuals within the sample. 

All diversity indices were calculated as log base 
2 values from the software BioDiversity Pro, Ver. 2 
(McAleece et al. 1997). Physico-chemical parameters 
such as water temperature, air temperature, pH, con- 
ductivity, salinity, total dissolved solid, and dissolved 



oxygen were analysed by Oakton waterproof Multi- 
parameter Tester PCS 35 and Oakton portable DO 
Meter. 

Structural diversity was calculated using statis- 
tical methods recommended by Heywood (2004) for 
floristic and taxonomic studies. We established the 
correlation of species composition and environmen- 
tal variables using Canonical Correspondence Analy- 
sis (CCA) with CANOCO Program for Windows 4.5 
package (Ter Braak and Smilauer 2002). Statistical 
significance of each variable was assessed using the 
Monte Carlo unrestricted permutation test involving 
999 permutations (Ter Braak 1990). The abbreviated 
names of species are given in the taxonomic table. The 
CCA plots represent overlap of species in relation to a 
given combination of environmental variables in the 
studied month. Arrows represent environmental vari- 
ables, with the maximal influence value for each vari- 
able located at the tip of the arrow (Ter Braak 1987). 
The Multiple Regression Analysis indicated a range of 
fluctuations in the algal community relative to the en- 
vironmental variables by construction of a correlation 
matrix between phytoplankton density and environ- 
mental variables and stepwise regression analysis with 
IBM SPSS 19.0, Bray Curtis cluster analysis with Bio- 
Diversity pro Ver. 2.0 program (McAleece et al. 1997) 
and comparative floristic studies with GRAPHS pro- 
gram (Novakovsky 2004). 

Results 

A total number of 44 species of phytoplank- 
ton of five taxonomic divisions, Charophyta (13), 
Cyanobacteria (12), Bacillariophyta (8), Chlorophy- 
ta (7), and Euglenozoa (4) were identified from the 
pool from March 2010 to February 2011 (Table 2). 
The monsoon season is represented by the maximum 
number of phytoplankton taxa (42), and the pre-mon- 
soon season with the lowest number (34). Among the 
44 phytoplankton taxa, charophytes were dominant 
throughout the study period. The percentage compo- 
sition of the studied phytoplankton (Fig. 2) reflects the 
uniform pattern of their representation (Charophyta 
> Cyanobacteria > Bacillariophyta > Chlorophyta > 
Euglenozoa) during the whole study period. The phy- 
toplankton density was highest in summer followed 
by the post-monsoon, monsoon and pre-monsoon 
densities (Table 3). The Shannon- Weaver diversity 
index (H') value (3.618) and the Pielou evenness (}') 
value (0.984) were found to be highest during the 
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Table 2. Diversity and ecology of algal species in plankton of the Shibpukur pool of the Burdwan district in March 2010 - February 201 1 



No 



Taxa 



Sub 



Oxy D Sal pH S Sap Aut Tro 



Cyanobacteria 

. Dolichospermum circinale (Rabenhorst ex Bornet & Flahault) P. Wacklin, L.Hoffmann p 

& J. Komarek 

2 Aphanocapsa grevillei (Hassall) Rabenhorst B,S 

3 Aphanocapsa litoralis (Hansgirg) Komarek & Anagnostidis P-B 

4 Aphanocapsa pulchra (Kutzing) Rabenhorst - 

5 Chroococcus turgid us (Kutzing) Nageli P-B,S 

6 Cylindrospermum bengalense Biswas - 

7 Merismopedia elegans A.Braun ex Kutzing P-B 

8 Merismopedia glauca (Ehrenberg) Kutzing P-B 

9 Merismopedia punctata Meyen P-B 

10 Phormidium formosum (Bory de Saint-Vincent ex Gomont) Anagnostidis et Komarek P-B,S 

11 Spirulina major Kutzing ex Gomont P,S 

1 2 Spirulina subsalsa Orstedt ex Gomont B 
Bacillariophyta 

13 Discostella glomerata (H.Bachmann) Houk et Klee - 

14 Eunotia sp. - 

1 5 Navicula trivialis Lange-Bertalot 

16 Nitzschia linearis (C.Agardh) W.Smith 

17 Nitzschia obtusa W.Smith 

1 8 Pinnularia major (Kutzing) Rabenhorst 

19 Surirella capronii Brebisson ex F.Kitton 

20 Ulnaria capitata (Ehrenberg) PCompere 
Euglenozoa 

21 Euglena viridis (O.F.Muller) Ehrenberg 

22 Phacus acuminatus Stokes 

23 Phacus longicauda (Ehrenberg) Dujardin 

24 Phacus pleuronectes (O.F.Muller) Nitzsch ex Dujardin 
Chlorophyta 

25 Chlamydomonas globosa J.W.Snow 

26 Chlorella vulgaris Beijerinck 

27 Coelastrum sphaericum Nageli 

28 Crucigenia tetrapedia (Kirchner) Kuntze 

29 Desmodesmus armatus (R.Chodat) Hegewald 

30 Monoraphidium contortum (Thuret) Komarkova-Legnerova 

31 Pediastrum duplex Meyen 
Charophyta 

32 Staurodesmus quiriferus (West & G.S.West) Teiling - 

33 Cosmarium angulatum (Perty) Rabenhorst - 

34 Cosmarium connatum Brebisson ex Ralfs - 

35 Cosmarium javanicum Nordstedt - 

36 Cosmarium quasillus PLundell - 

37 Cosmarium trilobulatum Reinsch - 

38 Desmidium baileyi (Ralfs) Nordstedt - 

39 Micrasterias foliacea Bailey ex Ralfs - 

40 Netrium oblongum (De Bary) Lutkemuller - 

41 Onychonema sp. - 

42 Pleurotaenium granuliferum (Joshua) Hirano - 

43 Spondylosium planum (Wolle) West et G.S. West P 

44 Staurastrum cingulum (West et G.S.West) G.M.Smith P 
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Note: Ecological types (Sub): B, benthic; P, planktic; P-B, planktic-benthic; pb, phycobiont; S, soil. Temperature (T): temp, temperate; 
eterm, eurythermic. Streaming and oxygenation (Oxy): st, standing water; st-str, standing- streaming. Saprobity (Watanabe et al. 1986) 
(D): es, eurysaprob; sp, saprophil. Halobity (Sal) (Hustedt 1938-1939): ph, polyhalobe; mh, mesohalobe; i, oligohalobious- indifferent; 
hi, oligohalobious-halophilous; hb, oligohalobious-halophobous. Acidity (pH) (Hustedt 1957): ind, indifferent; alf, alkaliphil; acf, acido- 
phil. Species-specific index of organic pollution (Sladecek 1986) (S). Saprobity (Sladecek 1986) (Sap): o, oligosaprob; o-b, oligo-beta- 
mesosaprob; b, beta-mesosaprob; b-o, beta-oligomesosaprob; b-a, beta-alfa-mesosaprob; a, alfa-mesosaprob; a-b, alfa-beta-mesosaprob; 
x, xenosaprob; x-b, xeno-beta-mesosaprob; o-a, oligo-alfa-mesosaprob; b-p, beta-polysaprob; i, eusaprob. Nitrogen uptake metabolism 
(Aut) (Van Dam et al. 1994): ate, nitrogen -autotrophic taxa, tolerating elevated concentrations of organically bound nitrogen. Trophic 
state (Tro) (Van Dam et al. 1994): me, meso-eutraphentic; e, eutraphentic. 
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Fig. 2. Dynamics of phytoplankton community composition as a percentage of taxonomic divisions 



post-monsoon season (Fig. 3). Simpsons dominance 
(D) index value (0.037) and Margalef s index of rich- 
ness (D ) (7.822) showed their highest values in the 
pre-monsoon season. 

The chemical variables slightly fluctuated dur- 
ing the studied period (Fig. 4) with the same ampli- 
tude of TDS, water and air temperature. The species 
numbers and abundance fluctuated in parallel, but op- 
posite to TDS and temperature, and correlated with 
the Shannon indices values (Fig. 5). 

An attempt was made to establish a correlation 
between chemical and biological variables (Table 4). 
Phytoplankton density showed a highly significant 
positive correlation with pH and salinity (p< 0.01), a 
significant positive correlation with air temperature, 
water temperature, dissolved oxygen and total sus- 



pended solids (p < 0.05), and a negative correlation 
with electrical conductivity (Table 5). In the regression 
analysis (Table 6) a mathematical model to ascertain 
the relationship between variables was proposed. The 
following model is used for the ordinary least square 
method regression: 

PHYTOPLANKTON DENSITY = a + (31 AIR 
TEMPERATURE + (32 WATER TEMPERATURE + (33 
pH + (34 SALINITY + (35 DISSOLVED OXYGEN + (36 
ELECTRICAL CONDUCTIVITY + (37 TOTAL DIS- 
SOLVED SOLIDS + e 

The response variable is phytoplankton diver- 
sity and the explanatory variables include air temper- 
ature, water temperature, pH, salinity, dissolved oxy- 
gen, electrical conductivity, and total dissolved solids, 
while e denotes the random disturbance term for the 
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Table 3. Diversity and abundance (cells ] 


oex litre) of algal 


species in 


plankton of the Shibpukur pool 


of the Burdwan district 


in March 


2010 


- February 2011. The abbreviation 


names are given for each species 


in the column co< 


ie. In the upper row as bold numbers are 


given 


the months of monitoring between March, 2010 as 


l,and 


i February, 2011 as 


12 
















No 


Taxa 


Code 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 




Cyanobacteria 




























1 


Dolichospermum circinale 


AnCir 


11 


11 


22 


2 


2 


5 


6 





11 











2 


Aphanocapsa grevillei 


AphaGr 


23 


22 


45 


12 


12 


10 


12 


8 


8 


8 


9 


10 


3 


Aphanocapsa litoralis 


AphaLi 


21 


22 


44 


11 


10 


12 


17 


11 


11 


20 


22 


22 


4 


Aphanocapsa pulchra 


AphaPu 


22 


33 


44 


11 





14 


17 


11 


11 


22 


22 


22 


5 


Chroococcus turgidus 


ChroTu 


22 


33 


55 


11 


12 


12 


14 


22 


22 


22 


22 


20 


6 


Cylindrospermum bengalense 


CylBen 


3 


5 


11 


2 


2 





1 


1 


3 


6 


6 





7 


Merismopedia elegans 


MerEle 


32 


32 


64 


16 





8 








8 


8 


16 


16 


8 


Merismopedia glauca 


MerGIa 


16 


16 


32 





16 


16 


8 


8 


8 


8 








9 


Merismopedia punctata 


MerPun 


16 


16 


32 








16 


8 


8 


8 


8 








10 


Phormidium formosum 


OscFor 


22 


12 


33 








1 


11 


11 


11 


13 


14 


14 


11 


Spirulina major 


SpiMaj 


11 


13 


33 


9 


11 


10 


15 


11 


12 





3 


9 


12 


Spirulina subsalsa 
Bacillariophyta 


SpiSub 


22 


22 


43 





2 


11 


13 


13 


14 


15 


11 


11 


13 


Discostella glome rata 


CycGlo 


2 


2 


2 


5 





3 


5 


5 


6 


6 


7 


15 


14 


Eunotia sp. 


Eunot 


12 


10 


12 


3 


3 


1 


6 


7 


8 


9 


9 


15 


15 


Navicula trivia lis 


NavLan 


11 


10 





2 


2 


1 


8 


8 


9 


11 


15 


15 


16 


Nitzschia linearis 


NitLin 





1 


1 


4 


3 


2 


1 


5 


7 


8 


8 


18 


17 


Nitzschia obtusa 


NitObt 


4 





2 


11 


3 


3 





5 


5 


6 


6 


22 


18 


Pinnularia major 


PinMaj 


3 


3 





1 





2 


4 


4 


1 


4 


7 


12 


19 


Surirella capronii 


SurCap 


3 


3 





1 


1 


1 


4 


4 


5 


4 


5 


11 


20 


Synedra capitata 
Euglenozoa 


SynCap 


3 


3 


4 


1 


1 


3 


5 


5 


8 





12 


12 


21 


Euglena viridis 


EugVir 


2 


2 


6 


7 


9 





3 


3 


2 


5 


5 


6 


22 


Phacus acuminatus 


PhaAcu 


3 


3 


3 


5 


14 


5 


6 


6 


3 


3 


5 





23 


Phacus longicauda 


PhaLon 


3 


3 





1 


8 


16 


5 


5 


5 


7 


11 


11 


24 


Phacus pleuronectes 
Chlorophyta 


PhaPle 


5 





11 


11 


16 


17 


12 


4 


5 


6 


7 


7 


25 


Chlamydomonas globosa 


ChlaGI 





12 


11 


10 


3 


4 


12 


14 


15 
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Fig. 3. Variation of Biological Indices during the studied period 
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Fig. 4. Environmental variables, algal species richness and density fluctuation, March 2010-February 2011 
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Fig. 5. Fluctuation of the Shannon index and Species Richness, March 2010-February 201 1 

Table 4. Details of physico-chemical parameters of the Shibpukur pool between March, 2010 and February, 2011. In the first row are 
given the numbers of the monitoring months between March, 2010 as 1, and February, 2011 as 12 



Parameters 






1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


Air temperature (AT) 


[•C] 




32.1 


34.2 


36.3 


36.5 


35.5 


34.3 


31.2 


32.1 


32.0 


30.1 


31.5 


32.3 


Water temperature (WT) 


[°C] 




32.1 


32.2 


34.4 


34.4 


33.1 


33.1 


29.2 


29.0 


30.1 


29.0 


30.0 


30.3 


pH 






8.3 


8.5 


8.5 


8.3 


8.2 


7.9 


8.0 


8.2 


8.0 


8.2 


8.4 


8.3 


Salinity (Sal.) 


[g i- 1 ] 




2.3 


2.4 


2.6 


1.6 


1.6 


1.2 


1.2 


1.2 


1.3 


1.3 


1.5 


1.6 


Dissolved oxygen (DO) 


[mgl 


1 ] 


3.53 


3.46 


4.05 


2.66 


3.41 


3.47 


4.25 


4.22 


4.44 


4.15 


4.10 


4.07 


Electrical conductivity (EC) 


[uS cm- 1 ] 


7.8 


7.9 


7.5 


8.5 


8.5 


9.1 


8.2 


8.5 


9.0 


9.0 


9.2 


9.1 


Total Dissolved Solids (TDS) 


[mgl 


- 1 ] 


4.0 


4.4 


5.2 


5.2 


4.6 


4.6 


3.6 


3.1 


3.1 


2.9 


2.5 


2.6 



Table 6. Multiple regression model enumeration 



Explanatory variables 



Standardized coefficients (Beta) 



Air temperature (AT) 
Water temperature (WT) 

PH 

Salinity (Sal.) 

Dissolved oxygen (DO) 

Electrical conductivity (EC) 

Total Dissolved Solids (TDS) 

R Square Value 

F Value 



0.142(0.610) 

-0.694*** (-2.088) 

0.263(1.972) 

1.231* (6.672) 

0.554* (4.989) 

0.670** (3.558) 

0.835*** (2.457) 

0.986 

39.151 



Figures in parenthesis indicate t values; *, ** and *** indicate level of significance at 1, 5 and 10 percent respectively 
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Table 5. Correlation matrix among the physico-chemical parameters and phytoplankton density (Dens., in cells dm 3 ) in the Shibpukur 
pool of the Burdwan district between March, 2010 and February, 2011. (AT - Air temperature (°C), WT - Water temperature (°C), Sal. 
- Salinity (g l 1 ), DO - Dissolved oxygen (mg l 1 ), EC - Electrical conductivity (|^S cm 1 ), TDS - Total Dissolved Solids (mg l 1 )) 







AT 


WT 


PH 


Sal. 


DO 


EC 


Dens. 


TDS 




Pearson Correlat. 


1 
















AT 


Sig. (2-tailed) 
N 


12 


1 
















Pearson Correlat. 


0.933** 
















WT 


Sig. (2-tailed) 
N 


0.000 
12 


12 
















Pearson Correlat. 


0.293 


0.300 


1 












PH 


Sig. (2-tailed) 
N 


0.355 
12 


0.343 
12 


12 














Pearson Correlat. 


0.371 


0.502 


0.739** 


1 










Sal. 


Sig. (2-tailed) 


0.235 


0.097 


0.006 














N 


12 


12 


12 


12 












Pearson Correlat. 


-0.693* 


-0.758** 


-0.168 


-0.206 


1 








DO 


Sig. (2-tailed) 


0.012 


0.004 


0.603 


0.521 












N 


12 


12 


12 


12 


12 










Pearson Correlat. 


-0.414 


-0.446 


-0.450 


-0.751** 


0.207 


1 






EC 


Sig. (2-tailed) 


0.181 


0.146 


0.142 


0.005 


0.518 










N 


12 


12 


12 


12 


12 


12 








Pearson Correlat. 


0.105 


0.175 


0.720** 


0.863** 


0.226 


-0.518 


1 




Dens. 


Sig. (2-tailed) 


0.745 


0.587 


0.008 


0.000 


0.479 


0.084 








N 


12 


12 


12 


12 


12 


12 


12 






Pearson Correlat. 


0.885** 


0.909** 


0.130 


0.420 


-0.715** 


-0.612* 


0.074 


1 


TDS 


Sig. (2-tailed) 


0.000 


0.000 


0.686 


0.174 


0.009 


0.034 


0.819 






N 


12 


12 


12 


12 


12 


12 


12 


12 



** 



- Correlation is significant at the 0.01 level (2-tailed), * - Correlation is significant at the 0.05 level (2-tailed) 



proposed model. Standardized (3 coefficients for each 
explanatory variable with their significance levels are 
described in Table 6. The R 2 value (0.986) has estab- 
lished the applicability of the proposed model as a fit 
one. 

The CCA biplot (Fig. 6) was created for phyto- 
plankton species along with environmental variables. 
It showed specific zones of separation depending on 
special preferences of phytoplankton taxa with envi- 
ronmental variables. Circle 1 marks the species which 
prefer low density assemblages. These include Cos- 
marium quasillus, Desmidium baileyi, and Phormidium 
formosum. Circle 2 denotes species which are indica- 
tors of acidification, including Merismopedia elegans 
and Coelastrum sphaericum. Species of circle 3 are im- 



pacted by high values of dissolved oxygen and can be 
used as indicators of oxygenation. The species placed 
in this circle include Merismopedia punctata and Mer- 
ismopedia glauca. The triplot of CCA analysis (Fig. 7) 
shows phytoplankton species' abbreviated names and 
environmental variables. Arrows point to the month 
in which the marked variables are the highest. For ex- 
ample: DO is highest in Sept, Oct, Nov, Dec, and May 
(9, 10, 11, 12, and 5). 

Bio-indication of species responses to the 
changes in environmental variables shows that the 
pool community preferred low alkaline, low mineral- 
ized and low organic polluted water (Table 2, Fig. 8). 
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Fig. 6. CCA biplot with species abbreviated names and environmental variables. In circle 1 are marked species that are impacted by spe- 
cies richness. They preferred respectable low species numbers. In circle 2 are placed species that are impacted by high pH and therefore 
can be used as indicators of acidification. In circle 3 are placed species that are impacted by high DO and can be used as indicators of 
oxygenation. Species events are expressed as circles. Environmental factors are shown as arrows with the origin at their average values 
and extending towards higher values. pH - activity of the hydrogen ion, No sp. (number of species) - algal species richness, EC - electri- 
cal conductivity, WT - water temperature, TDS - total dissolved solids, DO - dissolved oxygen 
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Fig. 7. CCA triplot with species' abbreviated names and month of sampling (as a number of months). Arrows point to the month in 
which the marked variable is highest. For example: DO is highest in Sept, Oct, Nov, Dec, and May (9, 10, 11, 12, and 5). DO is higher 
when water temperature is lower because directions in which arrows of DO and WT are pointing are opposing. Species richness slightly 
correlated with increasing of WT, Density and DO, but with lowest EC, WT and TDS. No extremely high Species richness during the 
studied period, communities look like a continuum 



18 
16 
14 

m 
% 10 

m 8 
z 

4 
2 








































P P-B 

Substrate group 



B 



b 



10 




£ 6 

o 












LJ *» 

z 










2 











^_ 1 



acf ind 

Acidity group 



alf 




hb i hi mh 



Salinity group 



ph 



16 
14 



12 

a> 
5 10 

$■ 6 



c 
c 



d 



6 
4 
2 




0^" 



# 



• 



• 



^ 



# 






Class of Water Quality 



Fig. 8. Phytoplankton species composition as an 
indication of Substrate preferences (a), Acidity (b), 
Salinity (c), and Water Quality Class (d) 
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Discussion 

This study of the tropical climate community 
shows two-period seasonality as was found in the 
Eastern Mediterranean freshwater ecosystems (Bari- 
nova et al. 2010; Barinova 201 1), in which rainy season 
communities were more diverse and abundant (Laskar 
and Gupta 2009) in contrast to the diverse Central Eu- 
ropean small lakes (Dembowska 2011; Jaworska and 
Zdanowski 2011). In distinction, the moderate cli- 
matic zone phytoplankton communities show three 
periods of biomass production, which are related with 
the ice forming process (Dembowska 2011). Fluctua- 
tion of phytoplankton biomass and species richness 
in tropical lakes has a considerably smaller amplitude 
than in boreal lakes (Lewis 1990). Therefore, fluctua- 
tion of diversity in the Shibpukur pool with 44 species 
is similar to that in such pools as the Santragachi pool 
of North-East India in the Ganga River valley with 29 
species (Ghosh et al. 2012) and Chatla Lake with 34 
species (Laskar and Gupta 2009) and can be related to 
monsoon influence. 

The Shannon-Weaver index was used to exam- 
ine diversity in a categorical way which actually estab- 
lishes the information entropy of the species distribu- 
tion. This index considers the number of species and 
evenness of their distribution: an increase in evenness 
corresponds to a greater Shannon-Weaver value than 
there should be, which was supported by our results. 
This index also reveals the trophic status and pollution 
load of the water body. Our calculation of the Shannon 
index is similar to that for Chatla Lake (Laskar and 
Gupta 2009) and Santragachi Lake (Ghosh et al. 2012) 
in the Ganga River valley. The water body under con- 
sideration showed almost clean (Wilhm and Dorris 
1968) to slight polluted status (Staub et al. 1970). Simp- 
sons dominance index corresponds to the number of 
species present, along with the relative abundance of 
each species. Indeed, Simpsons index captures the 
variability of species abundance distribution. There is 
an opposite relationship between species dominance 
and diversity: when dominance increases the diversity 
decreases. As can be seen from our calculation, when 
there was maximum diversity value (3.618) in January, 
the dominance was the lowest (0.027). 

The Canonical Correspondence Analysis (CCA) 
biplot (Fig. 6) shows that dissolved oxygen (DO) in 
the pool is higher when the water temperature (WT) 
is lower: the arrows of DO and WT point in opposite 
directions. Species richness is weakly correlated with 



an increase in WT, density (Dens.), and DO, but is in- 
versely correlated with electrical conductivity (EC), and 
total dissolved solids (TDS). As was studied in Chatla 
and Santragachi Lakes in similar climatic condition, 
EC, WT and TDS are also major variables that impact- 
ed phytoplankton species diversity (Laskar and Gupta 
2009; Ghosh et al. 2012). While in Chatla Lake 34 algal 
species have been found during the year-round moni- 
toring study, we found 44 species in the Shibpukur pool 
phytoplankton communities during the same period. 
Both lakes do not show well-expressed species rich- 
ness peaks during the studied period, and the phyto- 
plankton communities are stable in respect to species 
content. Our CCA analysis thus reveals indicators for 
species richness, acidification and oxygenation. These 
relationships predict future changes in the ecological 
status of this type of water bodies. 

The bio-indication method, implemented for 
the first time in a study of phytoplankton in Indian 
pools, shows three types of algal habitations: 1) plank- 
tonic, 2) planktonic-benthic and 3) benthic-periphy- 
tonic. In a small shallow pool, the attached and ben- 
thic forms mix with planktonic forms as a result of 
monsoon water disturbance. 

The comparative floristic study, performed for 
the first time in Indian pools, allows clusters of phy- 
toplankton taxa to be distinguished in respect to their 
taxonomic similarities and phylogenetic affinities. A 
similarity tree of floristic composition was constructed 
for the phytoplankton communities of a protected pool 
in the Burdwan district (Fig. 9), showing two clusters 
at the similarity level of 40%. The first cluster includes 
the phytoplankton assemblages collected in May, June, 
July, October, January, February and the second cluster 
comprises the assemblages for the rest of the year. The 
dendrite of floristic (Fig. 10) and Bray-Curtis analysis 
(Fig. 11) using the complete matrix of monthly phyto- 
plankton communities represent similarity patterns of 
the two clusters of phytoplankton assemblages. 

Conclusion 

This paper examines the detailed information on 
phytoplankton composition and environmental factors 
and the temporal pattern of their representation in a 
stagnant water body that is typical of pools in North- 
Eastern India, which are influenced by the monsoon cli- 
mate. Freshwater ecosystems are subject to various en- 
vironmental impacts which cause changes in their phy- 
toplankton composition (Cetin 2000) in particular af- 
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Fig. 9. Similarity tree of phytoplankton floristic composition (50% similarity level) reflects two different complexes. Here the bold 
number indicates the month starting from March 2011 as 1 to February 2011 as 12 




Fig. 10. Comparative species list statistics as dendrite of similarity (50% similarity level) of monthly communities. The bold number 
indicates the month starting from March 2011 as 1 to February 2011 as 12 
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Fig. 11. Bray-Curtis cluster analysis (complete link) for establishing monthly community similarity of phytoplankton 



fecting tropical lakes' phytoplankton dynamics (Melack 
1979). The Ganga River valley pool communities exhib- 
it maximal diversity of phytoplankton taxa in the mon- 
soon season. The comparison of 44 phytoplankton spe- 
cies from the Shibpukur pool, an undisturbed stagnant 
water body in a protected area, reveals a similar pattern 
of species representation throughout the studied peri- 
od. The successive communities form a continuum cor- 
responding to the Constancy (C) category of (Colwell 
1974). The continuum probably reflects low amplitude 
environmental fluctuations, stagnant water conditions 
and the absence of any anthropogenic or other distur- 
bances except the monsoons. Because the Shibpukur 
pool phytoplankton communities are similar in vari- 
ation of algal abundances and environmental variables 
to many other pools of the Ganga River valley (Laskar 
and Gupta 2009; Ghosh et al. 2012), but studied in more 



detail with regard to species richness, abundance, bio- 
mass, biological indices, CCA, comparison of commu- 
nities, and bio-indication, it can be used as a model of 
aquatic community dynamics under seasonal fluctua- 
tion in the monsoon climate, applicable for monitoring 
of water bodies in the West Bengal Province. 
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